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Abstract. We present the Spectral Energy Distribution of HR 4049 based on literature data and new continuum 
measurements at 850 fim. The SED shows variable absorption in the UV, and a large IR excess, both caused by 
circumstellar dust. The shape of the IR excess from 1 fim all the way down to 850 fim can be nearly perfectly 
fitted with a single blackbody function at T « 1150 K or alternatively with a sum of blackbodies in a narrow 
temperature range. The energy emitted in this IR continuum radiation is about one-third of the stellar luminosity. 
We show that this blackbody radiation must be due to the presence of a circumbinary disk with a large height. 
This disk must also be gas-rich, in agreement with the observations of molecular bands in the ISO-SWS spectrum. 
We present two possible scenario's for explaining the shape and the intensity of the IR excess. The first scenario 
involves large grains (a > 1 mm) that each radiate like a blackbody. The second scenario argues that the blackbody 
radiation is due to a very optically thick circumbinary disk. We investigate if such a disk would indeed produce 
blackbody radiation by presenting results from radiative transfer calculations. We further quantify the properties 
of such a disk and its stability in the framework of (hydro)dynamics, grain settling, radiation pressure and grain 
drift. The virtues and shortcomings of both models for the origin of the IR blackbody are discussed by contrasting 
them with other observations and assessing them in the framework of (binary) (post-)AGB evolution. 
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1. Introduction 

When low- to intermediate mass stars leave the 
Asymptotic Giant Branch (AGB), they spend a relatively 
short time (10'^ to 10'* years) on their way to become 
a planetary nebula; this stage of stellar evolution is the 
post-AGB phase. During this phase, the stellar interior 
is gradually more exposed, so that the effective temper- 
atures increase from typical AGB values (~ 3000 K) to 
those commonly found in planetary nebulae (~ 30 000 
K). At the same time, the luminosity of these stars re- 
mains more or less constant and is generally in the range 
10^ - 5x lO'^ L0. As the AGB phase is characterized by 
substantial mass-loss, it is expected that post-AGB stars 
are surrounded by circumstellar matter. It is this ma- 
terial that might later on become visible as a planetary 



nebula. 

HR4049 is considered to be the prototype of a group of 
post-AGB stars sharing similar - but rather peculiar - 
characteristics (see van Winckel et al., 1995). They are 
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all members of a (close) binary system, often with highly 
eccentric orbits. Their photospheres show a high degree 
of metal-depletion, while volatile elements as CNO, S and 
Zn show solar-like abundances. Most of them exhibit a 
deficiency in the UV and optical, while they do have a - 
sometimes huge - IR excess, pointing to the presence of 
circumstellar dust which thermally re-emits the radiation 
it absorbs in the UV. The key to explaining all of these 
peculiarities is probably the binary nature of these objects 
which allows the formation of a massive, stable and long- 
lived circumbinary disk. The formation of dust in these 
disks followed by re-accretion of the depleted gas onto the 
stellar ph otosphere expl a ins th e photospheric abundance 
patterns ( Waters et al. , 1992| ). Moreover, the interplay 
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of a massive circumbinary disk with the central binary 
might be important for the dynamics of the binary sys- 
tem ( |Mellem4 |1997| ; |Icke| , |198^ ). The lead role of the dust 
disks in these systems can therefore hardly be overesti- 
mated. 

HR4049 was first suggested to be an object in the transi- 
tion phase from the AGB to the Planetary Nebula stage 
by Lamers et al. (1986) who discovered both a large IR 
excess and severe UV deficiency indicating the presence of 
circumstellar dust. From a variety of observational results, 
it is clear that the properties of this circumstellar dust are 
unusual. 



Waelkens et al. (1991) showed that photometric variations 
are correlated with radial velocity variations : the star is 
faintest and reddest when the star is at inferior conjunc- 
tion. The most plausible explanation for this phenomenon 
is that the dust which is responsible for the extinction is 
not distributed in a spherical shell, but concentrated in a 
disk surrounding the entire binary system - a circumbi- 
nary disk - observed at a rather high inclination : in such a 
geometry the amount of dust in the line of sight is largest 
at inferior conjunction which explains both the photome- 
try and color variations. 

The amplitude of the circumstellar extinction is 
wavelength-dependent, and scales linearly with A^^ 



(Waelkens et al., 1991; Waters et al., 1989; Buss et al 



1989). The extinction remains linear in the UV and shows 
no sign of the 2175 A bump (Buss et al., 1989] ). A lin- 



ear spectral extinction can be understood in terms of Mie 
scattering when the size of the dust grains is small com- 
pared to the wavelength. This yields an upper limit of 300 
A for the size of the dust grains in the line-of-sight towards 
HR4049. |Joshi et al.| ( [1987D and [Johnson et aL| ( |1999D pre- 
sented spcctropolarimetry in the UV and the optical. The 
data show a high degree of polarization in the UV, with 
a change of 90° in the position angle around 2000 A- the 
signature of scattering material in a bipolar distribution. 
The p olarization in the opt ical seems to vary with orbital 
phase ( Johnson et al , |l999| ), suggesting that the polariza- 
tion in the optical is due to scattering in the circumbinary 
disk. In the UV, the polarization is then caused by scatter- 
ing in the bipolar lobes, which should contain a population 
of small grains (a < 0.05 /im). 



Waters et al.| ( 1989 ) studied the properties of the cir- 
cumstellar dust around HR4049. They detected emission 
features attributed to Polycyclic Aromatic Hyd rocarbons 
(PAHs 
et al 



the ISO/SWS spectrum of this object ( Bcintema 



around an 0-rich AGB star. The dust surrounding HR 
4049 therefore seems to be C-rich, while the gas is more 
0-rich in nature. 

Another intriguing property of the ISO/SWS spectrum is 
that the entire IR continuum can be matched nearly per- 
fectly with a blackbody at a single temperature of 1200 
K. An extrapolation of this blackbody also matches the 
60 /zm IRAS point. 

In this paper we (re-)investigate the circumstellar envi- 
ronment of HR 4049. In Sect. ^ we present the Spectral 
Energy Distribution (SED) of HR 4049 from the UV to 
submm wavelengths, showing that the dust continuum re- 
mains a blackbody to the longest wavelengths observed 
so far. In Sect. § we show that the SED can be used to 
derive strong constraints on the circumstellar dust distri- 
bution. We conclude that there are only two possible dust 
distributions - both having a disk geometry - that can 
meet these observational constraints; these are described 
in Sect. ^ and ||. The properties of the circumstellar dust 
disk within the framework of these two possible scenarios 
are addressed in Sect. ^. In Sect. we compare the two 
possibilities in the framework of disk dynamics and stel- 
lar evolution, and assess the plausibility of both models 
in the framework of other observations presented in the 
literature. 

2. The Spectral Energy Distribution of HR4049 

To compose the Spectral Energy Distribution (SED) of 
HR 4049, we used UV, optical and near-IR observations 
found in the literature, data obtained with the Short 



Wavelength Spectrometer (SWS, de Graauw et alj, 1996) 



[Kessle: 



on board the Infrared Space Observatory (ISO 
[st al. , 1996) and new data obtained with the James Clerk 
MaxweU Telescope (JCMT). 

The UV part of the SED is based on the many mea- 
surements of HR 4049 by the International Ultraviole t 
Explorer (lUE) (see e.g. Monicr fc Parthasarathy , 1999). 
Ground-based optical and near-IR photometry observa- 
tions were acquired over a long period; the observations 
we use in this paper are basically the same as presented in 
e.g. [van Winckej (|200l| ). In the IR, we use the ISO/SWS 
spectrum from 2.3 to 30 /im (see Paper I) and IRAS 
photometry points at 12, 25, 60 and 100 /im. Note 
however that the 100 /im point is an upper limit; this 
observation is probably contaminated by interstellar 
cirrus clouds. The 60 /im point is not contaminated 



] 996) shows these features in great detail, suggestive because the emission at this wavelength is clearly a point 



of a G-rich dust composition. This is further corroborated 



1989) 



by Geballe et al. (1989) who studied the profiles of the 3.43 
and 3.53 /im emission features recently identified with vi- 
brational modes of hydrogen-terminated crystalline facets 
of diamond (Guillois et al., 1999| ). 



In an accompanying paper (Paper I hereafter, in prepa- 
ration) we analyze the molecular bands found in the 
ISO/SWS spectrum. Although the dust features present in 
the spectrum can all be attributed to C-rich components, 
the molecular species are typical for the composition found 



source (Waters et al. 
In this paper we report upon very recent observations 
carried out with the JCMT using the Submillimetre 
Common-User Bolometer Array (SCUBA) to obtain 
a continuum measurement at 850 fim. HR 4049 was 
observed 3 times (31 Jan 2001, 9 Oct 2001 and 23 
Nov 2001). The data were reduced using Uranus as a 
calibrator; the atmospheric opacity at 850 /im is fitted 
with a straight line at epochs between the observations of 
photometric calibrators. The total integration time was 
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Fig. 1. The SED of HR 4049 at photometric maximum. The diamonds show lUE measurements, optical and near-IR 
photometry; these measurements are de-reddened using E{B — V) = 0.1. IRAS observations at 12, 25, 60 and 100 ^m 
are shown as crosses. Note that the 100 ^m point is only an upper limit. The continuum measurement at 850 /^m 
based on the new SCUBA observations is represented as the square. The ISO-SWS spectrum between 2.3 and 27 /Ltm 
is plotted as the solid black line. The thick grey line shows a Kurucz model for the star with parameters Tcs = 7500K, 
log 5 = 1 and [Fe/H] = —4.5. This Kurucz model is scaled to match the SED at photometric maximum at 580 nm. 
The inset shows the SED at photometric maximum (diamonds) and photometric minimum (+) at UV and optical 
wavelengths, revealing (variable) extinction from circumstellar dust. 



2700 seconds and yields a flux of 10.4 ± 2.2 mJy. We 
included this point in the SED. 

As HR 4049 shows considerable photometric variations 



from the UV to the near-IR (e.g. Waelkens et al., 1991) 



we composed a SED for both photometric maximum and 
for photometric minimum. 

The best estimate for the distance to HR4049 is 
given by the Hipparcos parallax of 0.150 ± 0.064 
mas which corresponds to a distance range be- 
tween 467 and 1162 pc. To study the properties of the 
circumstellar dust, we first need to correct the SED for 
interstellar extinction. This correction is somewhat un- 
certain, as the extinction towards HR 4049 contains both 
a (variable) circumstellar and an interstellar component. 
As the color variations are in phase with the photometric 
and radial velocity variations, the circumstellar reddening 
is lowest at photometric maximum. To minimize the 
effects of circumstellar extinction, we therefore used the 



SED at photometric maximum to estimate the interstellar 
reddening. 

We estimated the interstellar component of the E{B — V) 
by using the average interstellar extinction law given by 
Ravage fc Mathigj ( |1979D to de-redden the observed SED at 



photometric maximum and subsequently comparing this 
de-reddened SED with a Kurucz model. Several authors 
derived stellar parameters for HR 4049 using different 
methods (see Bakkcr ct al., 1998, for an overview); we 



used these parameters and compared the de-reddened 
SED therefore with a Kurucz model with Tes = 7500K, 
logg — 1 and [Fe/H] = —4.5. We found that adopting 
an interstellar reddening of £'is=0.1 yields the best fit in 
the optical between the SED at photometric maximum 
and the Kurucz model. This value for the interstellar 
reddening toward HR 4049 agrees well with the reddening 
determined for nearby stars which was found to be 
0.09 < Eis{B-V) < 0.17 ( |Buss et al.| , |l989[ ). We adopted 
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the same value for the Eis{B — V) at photometric 
minimum. Fig. Q shows the resulting de-reddened SEDs 
at photometric maximum and minimum and the Kurucz 
model. Note that the plot shows a clear UV of the 
dereddened spectrum compared to the Kurucz 
model with appropriate metallicity, a clear sign of 
anomalous circumstellar extinction which is also 



38 %. We can also use the SED to get an estimate of the 
optical depth : 



variable with orbital phase (Waelkens et al 



1995 



Monier &: Parthasarathy, 1999). In Fig. ^ we have 



subtracted the stellar model from the SED in order to 
show only the IR excess. 



3. Constraints on the dust distribution 

3.1. Energy considerations 

The SED shows a UV deficiency and an IR excess, both 
caused by circumstellar dust. The circumstellar dust ab- 
sorbs the stellar radiation in the UV producing the UV 
deficiency. The IR excess is caused by grains which are 
heated by stellar radiation and re-emit this energy ther- 
mally. It is important to realize (and often overlooked) 
that the IR excess is produced by all emitting dust grains 
in the beam, while the UV deficit is only due to absorp- 
tion by dust grains in the line- of- sight. The source of the 
UV deficit and the IR excess are therefore not necessarily 
the same dust particles. 

To get a feeling for the amount of absorption in the UV 
and the IR emission of the dust, we calculate the amount 
of energy absorbed in the UV and emitted in the IR by nu- 
merically integrating (parts of) the SEDs and the Kurucz 
model and express these energies in units of the stellar 
flux. All integrations were done using a five-point Newton- 
Cotes integration formula. 

We first calculated the stellar flux by integrating the 
Kurucz model between 145 nm and 850 fim. The energy 
radiated at shorter and longer wavelengths is only a neg- 
ligible fraction (SS 10"'*) of the total stellar flux. To cal- 
culate the amount of energy absorbed, we should inte- 
grate the SEDs without the IR excess at both photomet- 
ric maximum and minimum in order to find the fluxes 
^max ^j^^ ^min^ respectively. The IR excess dominates at 
wavelengths longer than 1.2 /im; the stellar flux at these 
wavelengths however still contributes about 10 % of the to- 
tal stellar luminosity. Fortunately the extinction becomes 
small at those wavelengths. We therefore assumed that 
the intrinsic SED at these wavelengths is equal to the 
Kurucz model. We then find 

1 ^ = 0.04 (1) 



1 



F, 



= 0.38 



(2) 



where the superscripts refer to photometric maximum 
and minimum respectively. At photometric maximum, the 
dust in the line-of-sight thus absorbs 4 % of the stellar ra- 
diation while at photometric minimum this amounts to 



r = -\n-^^ (3) 

-'Kurucz 

If the optical depth were grey, we can use Eqs. (|l|) and (|2|) 
to find that Tmax ~ 0.04 and Tmi„ 0.48. H o wever , the 
extinction in the UV is not grey ( |Buss et al.| , |1989[ ) and 
therefore the optical depth in the UV will be larger than 
these grey values. We used Eq. || to get an idea about the 
actual (wavelength-dependent) optical depth and found 
that the optical depth reaches maximum values of r ~ 1.5 
- 2 between 180 and 200 nm at photometric maximum, and 
even larger values (r ~ 2 - 2.5) at photometric minimum. 
Integrating over the IR excess yields 



Fir 
F, 



0.32; 



(4) 



so the energy radiated by the dust in the IR is about one 
third of the energy radiated by the star at all wavelengths! 
It therefore seems that the UV deficiency at minimum 
light and the IR excess both amount to about one third of 
the stellar flux. Looking only at the observations at min- 
imum light, one might therefore be tempted to conclude 
that the dust is distributed spherically around the star 
and that the line-of-sight we probe with the UV absorp- 
tion is representative. However, the large variation of the 
UV extinction with phase in fact tells us that the similar 
numbers are purely coincidental. We will get back to this 
in the discussion. 



3.2. Dust Temperatures 

In Fig. 1^ we have subtracted the stellar model from the 
SED in order to show only the IR dust emission. It shows 
a shape remarkably close to a blackbody over almost 3 
orders of magnitude in wavelength, from 1.2 to 850 /im. 
As we showed in Paper I, a Planck function at a single 
temperature of T ~ 1200 provides an excellent fit to the 
IR continuum up to 60 /im. In this Section, we will further 
quantify and constrain this dust temperature. 
The first question to answer is if a single temperature 
blackbody can really provide a good fit to the continuum 
in the statistical sense. We therefore decided to perform 
a minimization. As the ISO-SWS spectrum shows a 
wealth of features due PAHs, nano-diamonds and molec- 
ular bands (see Paper I), we preferred not to use the 
spectrum for this purposes and only restrict ourselves to 
the near-IR photometry observations, the IRAS measure- 
ments and the SCUBA point. However, when comparing 
the L and M band observations (3.8 and 4.7 /im) with the 
ISO spectrum, it is seen that the fluxes of the photometric 
points are slightly higher than the ISO-SWS continuum; 
the photometric points are likely influenced by the strong 
emission from the CO2 stretching band which is centered 
around 4.2 /tm. Also the 12 /im IRAS point is not repre- 
sentative for the continuum (see Fig. |l|); in this case, the 
culprit is probably the strong 11.6 PAH feature. Given 
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Fig. 2. The IR-excess of HR 4049. The Kurucz model is 
subtracted from the SED shown in Fig. |l|. Errors on the 
near-IR photometry points, the IRAS measurements and 
the SCUBA point are indicated, but are generally smaller 
than the symbol size. Note that the "photometric" points 
at 3.8 /im and 10 /im were derived from the ISO spec- 
trum (see text for more details). The prominent features 
in the ISO-SWS spectrum are mainly due to PAHs, nano- 
diamonds, and CO2 (see Paper I). The thick grey hue 
shows a single blackbody at a temperature of 1150 K. 
The inset shows that this single temperature blackbody 
provides a perfect fit to the entire continuum, including 
the 60 /im and 850 /im measurements. 

these considerations, we decided not to use these photo- 
metric points. Instead, we used the ISO-SWS spectrum 
to obtain "photometric" continuum points around these 
wavelengths. We replaced the L band photometric point 
with the mean flux in the ISO-SWS spectrum between 
3.7 and 3.9 /im; this range seems to be free of spectral 
features. The 12 /im IRAS measurement was replaced by 
the mean flux between 10 and 10.2 /im. It is somewhat 
uncertain whether this is really a good continuum point, 
as the PAH features might still contribute to the flux at 
these wavelengths. It is however a much better continuum 
point than the IRAS 12 /im one. 

We next scaled a single temperature blackbody to the K- 
band flux (2.2 /im) and calculated xl (reduced x^)- We 
found that a single temperature blackbody can indeed ac- 
curately reproduce the continuum (x^ < 1) provided the 
temperature of this blackbody is between 1135 and 1160 
K with a nominal best fit around 1150 K {xt ~ 0-3). Fig. || 
shows that this single temperature blackbody does indeed 
convincingly reproduce the entire IR-submm continuum. 
If the IR emission is really produced by dust grains, it is 
unlikely that all the dust is really emitting at a single tem- 
perature. Rather, a distribution of dust centered at a cen- 
tral temperature Tc with a small temperature range AT 
around it would likely produce a continuum that mimics 
a single temperature blackbody. In order to get a feeling 
for the allowable range of temperatures, we therefore per- 
formed a second set of x^ minimizations where we fitted 




Terrp;erat J 



Fig. 3. xt contours for a sum of equally weighted black- 
body functions. The horizontal axis is the central tem- 
perature in this sum, the vertical axis is the temperature 
range around this central temperature; the AT = case 
are single temperature blackbodies. The white contour is 
the xt — ^ level, indicating that such a sum of blackbodies 
is indistinguishable from a single temperature blackbody 
even if the temperature range is as broad as 400 K. 

the observations with a sum of equally weighted blackbod- 
ies with temperatures between Tc — and Tc -f . The 
X^ contour map is shown in Fig. ^. As can be seen from 
this Figure, the allowable temperature range for such a 
sum of blackbodies is 



and 



T„in > 880if 



< li2bK 



(5) 



(6) 



corresponding to a temperature range of about 450 K, 
centered around 1100 K. Dust grains with a grey opacity 
and with a radial power-law number density distribution 
N{R) oc i?^" would imply that we have to give weights 
to the blackbody functions at different temperatures, with 
W{T) oc T2". We also made a x^ fit for a a=2 and find 
temperature limits of 730 and 1285 K. 

3.3. Geometrical constraints 

The main goal of this paper is to investigate which kind 
of dust distribution could reproduce the blackbody shape 
and the luminosity of the IR excess. Waters et al. (1989) 



actually modeled the IR excess with a spherical dust shell 
composed of small dust grains. This model could repro- 
duce the IR excess down to 100 /im, showing that a fine- 
tuned distribution of small dust grains can result into 
a blackbody like SED. However, in order to achieve the 
blackbody slope at long wavelengths, the dust distribu- 
tion must be very extended since the contribution of low- 
temperature grains is needed. In the model by Waters 
|et al. ( 1989D the distribution extended out to 2500R* (cor- 
responding to 500AU) in order to account for the 60 /im 
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flux. Since there is now evidence that the spectrum con- 
tinues as a blackbody down to at least 850 iim, the re- 
quired distribution would have to extend to temperatures 
as low 3K, while the interstellar radiation field alone al- 
ready heats dust to considerably higher temperatures. We 
will therefore not investigate such a fine-tuned distribu- 
tion of small grains in a spherical shell any further. 
There are two remaining possibilities to explain the black- 
body shape of the IR-excess. The first possibility would 
be that the individual grains that are responsible for the 
continuum emission are radiating like a blackbody them- 
selves. In that case, the grains must be large compared to 
all wavelengths present in the SED, and the temperature 
distribution relatively narrow. The second possibility is 
that the individual grains are not blackbody emitters; in 
that case we propose that the blackbody radiation is pro- 
duced by an optically thick disk surface which radiates 
only in a narrow temperature range. In the remainder of 
the paper, we will investigate these two possibilities and 
their consequences. 

Assuming that the IR emission is produced by a disk with 
height H above the midplane (total thickness of the disk 
is 2H), the dimensions can be derived from the constraint 
that one third of the stellar luminosity should be absorbed 
by the disk. If the disk is infinitely optically thick to the 
stellar radiation (like a solid surface) and located at a dis- 
tance R from the star, then the disk surface is 2ttR ■ 2H; 
this surface must then cover a solid angle of ^ x AttR^, so 
we find for the height of the disk 



H _ 1 
^ ^ 3 



(7) 



This is actually a lower limit for the height. If the disk is 
not optically thick but absorbs on average only a fraction 
1 — e^'^ of the stellar radiation, the height would have to 
be larger by a factor 1/(1 — e""^). 



H 

'r 



1 1 



3 1 



(8) 



As the only assumption that went into Eqs. (^ and (^) is 
that the dust is in a disk, the limit on the height is a 
major constraint on any disk model. 

4. Narrow temperature distribution of blackbody 
grains 

We will now investigate the possibility that the blackbody 
shape of the continuum is due to a collection of "large" 
grains that each radiate like a blackbody. "Large" grains 
in this case means that the size of these grains must be 
larger than all wavelengths present in the SED; therefore 
the grains must be at least mm-sized. 

4.1. Required disk structure 

Using the equation of thermal equilibrium of dust grains, 
we can translate the limits imposed on the temperature 



from Eqs. (^) and (||) into limits on the distance of the 
dust grains from the central star : 



1 -R* 
2"r 



1/2 



(9) 



Using a stellar effective temperature of 7500 K and an 
effective radius i?* of 47 (Bakker et al., 1998| ), we find 



that the dust grains must be between 3.5 and 8 AU from 
the star. 

There are a few caveats on these distance limits. Since 
these numbers have been calculated assuming that each 
dust grains receives the unattenuated radiation from the 
star, the outer radius is an upper limit. If we take into 
account that a fraction of the stellar radiation gets ab- 
sorbed, it can be shown that the outer distance decreases 
for small optical depth with a factor of e~'^/^; assuming an 
optical depth of r = 0.5, the outer limit would therefore 
go down to 6.2 AU. According to Eq. (||) the same optical 
depth would require the disk height to reach 2H = 1.7R. 
In this way, the radial extent of the disk can be traded 
against its height. For very low optical depths, the disk 
must be extremely high. For modest optical depths, the 
disk height decreases towards the limit given in Eq. (^, 
but the radial extent must be smaller. 
We can also estimate the disk mass. The total absorbing 
surface is Aabs = 2ttR ■ 2H which must be equal to the 
cross section of N grains: 2ttR ■ 2H — Nua^. The total 
dust mass is then 



Afdust = 2.85 • 10-^Afc 



R 



H 



Pdust 



4AU 1.4AU 2.7gcm-3 1mm 



(10) 



It is clear that the major difficulty of the model with large 
grains in a narrow temperature range is how to realize the 
very particular arrangement needed; 

1. The disk needs to have a fine-tuned radial optical 
depth between one and two over all scale heights. 
Optical depths less than one will fail to absorb enough 
stellar radiation. Optical depths larger than two will 
produce a temperature range which is too wide. 

2. A disk is required which is more extended in the ver- 
tical then in the radial direction 

The scale height of such a disk could be provided either 
by a gas-rich disk with a hydrostatic distribution of gas in 
which the grains are suspended, or dynamically in a gas 
poor disk. We will discuss both scenarios. 

4.2. Dust grains in a gas-rich disk 

If the disk were gas-rich, the same hydrostatic pressure 
which is required to hold up the vertical height would 
also act in the radial direction, attempting to spread the 
disk in this direction. Since in a Keplerian disk the radial 
component of the gravitational force is balanced by the 
centrifugal force, stability of the disk would also require 
significant sub-Keplerian rotation in the outer parts of the 
disk in order to balance the large gas pressure gradient in 
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the disk. The resulting shear motion leads to transport 
of angular momentum which would again spread out the 
disk in radial direction. Another concern is grain settling. 
The settling time for grains is a gas-rich disk is given by 
eq. (p6|). If we scale the dust mass of 2.85 • 10~^Mq with 
a nominal gas-to-dust ratio of 100, we can estimate a gas 
mass in the disk of 2.85-lO-'*M0. Using 3.5 AU and 6.2 AU 
for the radial extent of the disk, we find an average surface 
density of 30 g/cm^. At 4 AU, this corresponds to a settling 
time of 150 yr for 1 mm dust grains. Grain settling would 
half the disk height in this time scale. This estimate is in- 
dependent of the grain size since the settling time depends 
on the ratio of surface density over grain radius which is 
constant for the above estimate. A large height could be 
kept for a longer time only if the gas-to-dust ratio in the 
disk would be much greater than 100. 

It is therefore highly questionable if a cloud of large 
dust grains supported by a hydrostatic gas-rich disk could 
reach a stable configuration in accordance with the obser- 
vational constraints. 



that the thickness is given by a few times the diameter of 
the largest bodies in the disk. To counteract the collapse of 
the disk, massive external stirring would be needed, which 
would spread the disk in vertical and radial direction, in 
contradiction to the observed narrow range of tempera- 
tures. 

An interesting analogy for this model is the star 
HD 98800 which was described by Zuckerman & Bccklin 



(1993). This star is a binary system composed of two pre- 
main-sequence K dwarfs. The dust disk, located around 
component B, absorbs about 20% of the radiation of this 
component. It is easy to see that the large scale height 
required cannot be due to hydrostatic equilibrium (see 
section ^.l| ) since the stars arc not luminous enough. 
Therefore the most likely explanation for this system is 
that it is indeed a debris system. However, the same 
arguments limiting the lifetime of such a system apply 
(Zuckerman fc Becklin, 1993). Possibly, the disk in this 



system has been produced by a recent destructive event, 
or is maintained by the interactions of the binary system. 



4.3. Dust grains in a gas-poor disk 

The other possibility is that the disk is gas-poor. It would 
then be a structure resembling the asteroid belt or the 
Kuiper Belt in the solar system, or the debris disks found 
around main-sequence stars (e.g. Backman & Paresce, 
1995 ; [Chen fc Jura , 200l| ). In such disks dust grains are 



produced by collisions between larger bodies. Since grains 
smaller than about 1mm would be blown away by radia- 
tion pressure, this would also lead to a natural explanation 
for the absence of small grains indicated by the black- 
body spectrum. In a gas-poor disk, the dust grains move 
on Keplerian orbits around the central star. The vertical 
structure of such a disk is then given by the distribution 
of orbital inclinations. One could imagine that the radial 
extent of such a disk is set by a sheparding planet, much 
like it is in Saturns rings. 

However, it is easy to show that such a disk cannot 
keep the required large inclination distribution. As we 
have discussed above, the disk must have an optical depth 
of order one in the radial direction. Since the disk is higher 
than its radial extent, this immediately implies that the 
vertical optical depth will be of order one or even larger. 
Therefore, a dust grain moving on an inclined orbit will 
experience at least two collisions during each orbit. This is 
simply due to the fact that both the collision cross section 
and the absorption cross section are set by the geometrical 
cross section of the grains. 

The collisions between grains will either destroy the 
grains if the velocity is larger than a few km/s (Tielens 
et al., 1994), or dissipate energy. The energy dissipation 
leads to a rapid flattening of the disk. Jeffreys (1947) 
showed for the case of Saturns rings that such a disk 
quickly collapses to a mono layer thickness. A more de- 
tailed study including the effects of gravitational stirring 



5. An optically thick disk: The wall model 

The second explanation for the blackbody spectrum of 
HR4049 is that the emission comes from an optically thick 
disk surface. 

The scenario which we have in mind is a highly opti- 
cally thick cylindrical (see Fig. [l^ ) circumbinary dust disk 
which is illuminated from inside by the central star. The 
high optical depth prevents radiation to penetrate deeply 
into the disk which is therefore cold and invisible against 
the bright inner rim. What we observe as the blackbody 
radiation is the part of the inner disk rim which is directly 
visible to the observer. 

The disk will have a very hot inner rim, which will be at 
a single temperature if the rim is concave (i.e. at constant 
distance from the star). A vertical or more convex inner 
rim will produce a spectrum that is a bit more multi- 
temperature, but this deviation may still lie within the 



temperature constraints derived in Sec. 3.2 



The optical depth of this disk must be high, at all 
wavelengths considered. In fact, the optical depth must 
be so high that the colder dust which is not directly illu- 
minated by the star does not contribute to the emergent 
spectrum. We will quantify this below. 

In the optically thick model, the constraints on the 
particle size and on the radial extent which we derived in 
Sect. ^ can be relaxed. The individual grains do not have 
to be large, and there is no limit on the radial extent. 

5.1. Geometrical constraints revisited 



In Sect. 3.3 we have shown that the height of an optically 



of bodies embedded in the disk (Cuzzi et al., 197£) showed 



thick disk must he H ^ in order to absorb sufficient 
radiation, which was derived from the ratio of IR and stel- 
lar luminosities. 

We envision the emitting surface to be optically thick 
and radiating only into a half-space (see Fig. |l^. The flux 
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seen by an observer will depend upon the viewing angle 
under which the system is seen. From a nearly pole-on di- 
rection, the projected surface would be small and the IR 
radiation received low. From a nearly edge-on direction, 
the projected surface is large and the observed flux high. 
The observed IR luminosity therefore roughly scales with 
sini. Since the H/R ratio depends on the emitted IR lu- 
minosity rather than the observed luminosity, there is an 
uncertainty in the real value of H/R, depending on the 
inclination at which we observe the system. Because of 
the significant variations in circumstellar extinction dur- 
ing the orbital motion of IIR4049, a larger (near edge-on) 
inclination seems more likely, meaning that the real H/R 
could be somewhat smaller than i. 

From the temperature of the dust one can easily com- 
pute the distance of the inner rim to the central star. 
This distance is a factor of two larger than the distance 
of black body grains in an optically thin environment be- 
cause the surface can cool only into a solid angle of tt 
instead of An in the optically thin case. If this rim is lo- 
cated at a distance R from the star, then the stellar flux 
it receives is: = {R^/R)'^aT^. The flux emitted by the 
wall is F„aii — <jT^. Therefore the temperature of the 
wall, if located at a distance R/\sT = {R^/Rf/'^T^. For 
T = 1200K, ^ 7500K and i?, 47i?o one obtains 
a radius of the inner surface of 8.5 AU. If one considers 
the fact that the wall is also irradiating itself, then one 
must take into account that the total irradiation flux of 
the wall will be about a factor of 1 + fi/47r higher than 
the direct stellar irradiation. By taking ^1/Att ~ 1/3, as 
derived above, one finds R = ^/IJi x 8.5 AU = 9.9 AU. 
It should be noted that the estimate of 1 + Vl/An as the 
enhancement factor for the irradiation is very rough, but 
suffices for now. For a detailed discussion, see DuUcmon^ 



isoor 



et al l (^OOll) 



5.2. Emission from an irradiated dust wall 

We now turn to the question of whether an irradiated wall 
will indeed produce an emission that is sufficiently close 
to a blackbody curve to fit the SED of HR4049. Since 2D 
radiative transfer calculations with full frequency depen- 
dence and at optical depths reaching a million are very 
cumbersome, we will tackle this problem in a two-step 
approach.: First we will use a 1-D radiative transfer cal- 
culation in order to see if the inner surface of the wall does 
indeed produce blackbody-like emission. Then we will use 
a 2-D diffusion calculation in order to see if the emission 
of the rest of the disk (surface and outer rim) can be suf- 
ficiently suppressed using high optical depth. 

For the 1-D calculation, we consider infinite half-space 
filled with dust: for x < there is vacuum, for a: > 
there is a constant density medium of dust grains. This 
half-space is illuminated (on the x < side) with a flux 
F* which is taken to be the stellar flux as observed at a 
radius of 8.5 AU (we ignore the effects of self-irradiation 
mentioned above). We run the calculation for a mixture of 
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Fig. 4. Temperature structure of the ID calculation. 



grains of two sizes: small (0.1 /im) grains , and la rge grains. 
We take the opacity of Drainc fc Let (1984) for small 
grains, and a simple grey opacity for our 'large' grains. 
The grey opacity guarantees that the optical depth will 
be large at all wavelengths, a very important requirement 
of the wall model. We include the small grains in order to 
see if the superheating of small grains close to the surface 
of the disk will cause significant deviations from the black- 
body shape. The requirement of large optical depth at all 
wavelengths could also be met with non-grey, small-grain 
opacities. In this case the limiting criterion would be that 
the optical depth at 850 /im must be large. The transfer is 
performed using a frequency-dependent discrete-ordinate 
method (i.e. ray-by-ray transfer at each frequency bin and 
for each angle). The transfer equation is: 



dl. 



^^^p.<[i?.(r.)-/^,.] 



dx 



(11) 



where kJ, and are the opacities of the small grains and 
the large grains respectively, and idem for the densities 
pi and p2 and the temperatures Ti and T2. We ignore 
scattering for simplicity. The symbol /i is the cosine of 
the angle of the radiation. At every point x in the slab, 
the temperature of the dust grains is determined by the 
radiative equilibrium equation: 



poc poo 

/ KvBu{Ti)dv = / Ki,Judv 
Jo Jo 

where the mean intensity is given by 



I,j,,ud^i 



(12) 



(13) 



Using Lambda Iteration with Ng acceleration we it- 
erate this set of equations to convergence and find the 
solutions shown in Fig. ^ 

One sees that for the large grains, the temperature re- 
mains almost constant as a function of optical depth. The 
small grains superheat on the disk surface. At higher opti- 
cal depth, the two temperatures become the same because 
at large optical depth both grain species become radia- 
tively coupled. By performing direct ray-tracing through 
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Fig. 5. The SED produced by the ID calculation. The 
solid curve is indistinguishable from a single temperature 
blackbody. 



this model atmosphere, we can obtain the predicted SEDs, 
which are shown in Fig. |5[ The emerging spectra are close 
to pure blackbody curves, in fact the fat curve is indistin- 
guishable from a blackbody. The superheated small grains 
produce a slight broadening of the SED towards short 
wavelength. We conclude that an irradiated high optical 
depth surface indeed radiates like a blackbody and could 
in principle be used to explain the SED of HR4049. 



5.3. DifFusion though the disk 

In order to address the issue of diffusion through the 
disk towards the surface and towards the outer edge, 
we need to perform a 2-D radiative transfer calculation. 
Unfortunately, realistic radiative transfer calculations in 
2-D are rather cumbersome and complex, in particular at 
the very high optical depths needed for our disk model. 
Under these circumstances, a reasonably reliable approxi- 
mation of the 2-D transfer problem can be made using the 
equations of radiative diffusion. 

Since we are interested in the qualitative behavior of 
the radiative diffusion through the disk, we keep the disk 
geometry highly simplified. We assume a constant density 
torus limited in radius by 9AU < R < 18AU and in Q (the 
angle from the pole) between 7r/2— 1/3 < Q < 7r/2 -1-1/3- 
consistent with the H/R =1/3 requirement. We compute 
first how the direct stellar radiation penetrates the inner 
few mean free paths of the inner edge: 



F;(i?,e) 



47ri?2 



exp ■ 



pn^dR 



(14) 



We now assume that the absorbed radiation is re-emitted 
as IR radiation. This diffuse IR radiation is treated as 
a separate radiation field. We do this in a frequency- 
averaged way, and write down the diffusion equation for 
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25.0 



Fig. 6. Temperature contours of the diffusive disk. The 
temperature at the inner boundary is 1200K, further con- 
tours are in steps of 50K. 



the frequency-averaged mean intensity J(i?, 6): 

1 a / sin e 9 J 
i?sineae vs^p^ae 



i?2 OR \3pKdR 



-s 

(15) 



where k is the Rosseland mean opacity. This equation de- 
scribes the radiative diffusion in polar coordinates, derived 
from the first two moment equations with Eddington clo- 
sure. We also implicitly assume (for simplicity) that the 
mean opacity does not change very much from location 
to location, since otherwise additional terms have to be 
included in the diffusion equation. 
The source term 5" is given by: 



4tt 



pK^F*diy 



(16) 



This source term describes the transition of energy from 
primary stellar radiation into reprocessed diffuse infrared 
radiation. Since all primary radiation is absorbed in the 
innermost few mean free pathlengths, this source term will 
be only non- negligible close to i? = Ri^. So for larger radii 
the above diffusion equation (|l5| ) effectively states that the 
divergence of the flux should be zero, i.e. that energy is 
conserved. The boundary conditions at all edges are taken 



to be those described in Rybicki & Lightman (1979). 

The 2-D diffusion equation can now be solved by writ- 
ing it into a matrix equation, and solving this using stan- 
dard techniques such as the biconjugate gradient method 
(e.g. Numerical Recipes 1992). Once we find J{R,Q), we 
can compute the dust temperatures from: 



-T(i?, 9)4 = J(i?, 9) 

TT 



s 



pKp(T) 



(17) 



where Kp{T) is the Planck mean opacity at the dust tem- 
perature. 

We run several models, with an equatorial Rosseland 
optical depth between r — 100 and r — 10^ and one for 
T — 10^. The resulting temperature profile for the t — 10^ 
case is shown in Fig. |. InFig.|7| we show the SEDs of these 
tori, as seen from an inclination of 45 degrees. We used a 
2-D ray-tracing utility to produce these SEDs from the 
temperature distribution resulting from the 2-D diffusion 
equations. 
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Fig. 7. Fits to the observed SED with diffusion calcula- 
tions for tori with different optical depth. The solid grey 
line shows the 1150K blackbody fit. The dashed lines are 
(from top to bottom) for optical depth 10^, 10"^, 10^, and 
10^. 

One sees that even for very high optical depth the spec- 
trum will exceed a pure blackbody curve at long wave- 
lengths. The xl for the r = 10^ fit is 2.85. For high opti- 
cal depths only a very tiny amount of radiation is leaking 
through the disk. But since the Rayleigh-Jeans slope of 
a blackbody curve does not contain much energy, even 
this tiny amount of leakage shows up at long wavelengths. 
These results are in agreement with the simple estimate 
that the flux at the surface and the outer edge of the 
torus is proportional to l/r, where r is the optical depth 
towards the inner edge. If we assume that the flux is pro- 
portional to T'^, then one can estimate that the temper- 
ature of the outer- and upper/lower surfaces of the torus 
is proportional to Tsurf oc t~^/^. This means that in order 
to suppress the temperature of these back-surfaces of the 
disk by a factor of 0.1 one must have an optical depth of 
T ~ 10^. A pure blackbody emission from this torus model 
therefore requires very high optical depth. 



5. Will the disk remain radially mixed, or will the radia- 
tion pressure of the star lead to a grain-gas separation? 

6. Will radiation pressure of the star be sufficient to blow 
away (parts of) the disk? 



As we already discussed in Sect. 4.3, a gas-poor disk 
cannot keep the required scale height because of frequent 
collisions. This argument is even much stronger for an op- 
tically thick disk. We will therefore assume the disk is 
gas-rich and in hydrostatic equilibrium, to which we can 
apply the standard theory of hydrostatic disks. We will 
adopt a cylindrical axisymmetric coordinate system with 
radial coordinate r and vertical coordinate z. We will often 
refer to scale heights in the following sections; we remind 
the reader that these scale heights correspond to the verti- 
cal coordinate starting from the equatorial plane. For the 
following estimates we will for simplicity assume that the 
disk has a constant surface density Sq and a constant scale 
height ho (i.e. a non-flaring disk, see below) in a distance 
range Ar — ri — ro from the star. We will further assume 
that inside the disk the dust-to-gas ratio is /d = 0.01 and 
that the vertical distribution of the gas is given by hydro- 
static equilibrium in an isothermal disk. We assume the 
dust particles all to have the same size a and a specific 
density pd- 

Armed with these properties, we may calculate the to- 
tal disk mass 



(18) 



When evaluating equations numerically in the follow- 
ing, we will always use a standard model of the disk which 
contains O.SM©, with tq = 10 AU and n = 20 AU. This 
corresponds to a surface density of 1414 gcm~'^. The disk 
mass of O.3M0 is very high and should be considered 
an upper limit for the possible disk mass since signifi- 
cantly higher masses will lead to gravitational instabilities 
if A/disk > ( Pringlc , 1981 ). We use this high mass in 

order to maximize the stability of the disk against grain 
settling or radiation pressure. 



6. Disk properties 

The disk we are proposing to be the source of the 1150 K 
blackbody radiation in IIR4049 certainly is unusual, be- 
cause of the high scale height (covering fraction) at the in- 
ner boundary and the high optical depth. In this Section, 
we will derive properties of the disk including mass esti- 
mates and stability considerations in order judge how re- 
alistic such a model can be. We will therefore investigate 
the following questions. 

1. Can we understand the large disk height {H/ R = 1/3)? 

2. Can we place limits on the grain size? 

3. What disk mass is required to make the disk suffi- 
ciently optically thick? 

4. How long will it take for the dust grains to settle to- 
wards the midplane? 



6.1. Vertical scale height from hydrostatic equilibrium 

We assume the disk to be gas-rich. Its height can therefore 
be calculated from the assumption of hydrostatic equilib- 
rium. The disk pressure scale height hp is given by the 
square root of the ratio of the sound speed an d the local 
Kepler frequenc y and can be written as (e.g. Chiang fc 
|Goldreich| , |l997| ) 



hp / kTr 

G{Mi + M2)fimp 



(19) 



where Mi and M2 are the masses of the two stars in the 
IIR4049 system, /i is the mean molecular weight and nip 
is the proton mass. It should be noticed that this formula 
holds for an annulus of a thin disk, and may be inaccurate 
for the inner edge of a thick disk. For r = 8.6AU, Mi ~ 
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0.56 Mq and M2 ~ Mi, and T = 1150K this amounts to 



^ = 0.19 



(20) 



If we compare this computed scale height with the ob- 
served required scale height (see Sect. O), which was 
H/r = 1/3, and if we take into account that the surface 
height (the height at which t — 1) is usually a few times 
the pressure scale height we arrive at 



H H 
— = 0.19— 

r hp 



(21) 



where H/hp is roughly in bet ween 1 and 4 for a Gaussia n 
vertical density distribution ( Chiang fc Goldreich , 1997 ) 

Therefore, a hydrostatic disk can reach the required 
height without any difficulty. In fact, a mechanism to trun- 
cate the disk somewhere between one and two pressure 
scale heights seems to be required in order to match the 
observational constrains. 

In the passive disks around young stars, the surface 
irradiation of thin circumstellar disks leads to an increase 
of the ratio hp/r with distance from the star. Such disks 
are called flaring disks, and have been successfully used to 



explain the flat SEDs of young stars with disks (Kcnyon 



& Hart |mann| , |l987| ; [Chiang fc Goldreich| , |l997D . Could the 
disk around HR4049 also be flaring? The blackbody shape 
of the SED is clear evidence that this is not the case, be- 
cause flaring disks generally have very broad, multicom- 
ponent SEDs. The physical criterion for flaring to occure 
is that the disk midplane temperature must decrease more 
slowly than r~^, so that hp/r increases with distance (see 
eq. (p9|)). In the wall model, the disk temperature is not 
determined by the radiation reaching the disk surface. 
Instead, radiation is only impacting on the inner rim and 
transported outwards through the disk by diffusion. The 
temperature will therefore decrease quickly with distance 
from the star, and hp/r will decrease, preventing any flar- 
ing. The shape of such a diffuse disk inner rim has recently 
been discussed in detail in connection with the near IR 



emission of Herbig AeBe stars (Dullemond et al, 2001) 



6.2. Radial optical depth 

The density structure p{z) in the disk assuming hydro- 
static equilibrium is given by (Chiang & Goldreich, 1997) 



27r/io 



exp 



2hl 



(22) 



If dust and gas are perfectly mixed, the number density 
of dust grains is given by 



nd(z) 



(23) 



where fd is the dust-to-gas ratio and is the mass of a 
single dust grain. The dust opacity is then given by 



Cd{z) = ndQ{a, A)7ra^ 



Sfdp{z)Qia,X) 
4apd 



(24) 




Fig. 8. The optical depth for short wavelength radiation 
along radial rays through the torus as a function of height. 
The different curves are for different grains sizes with fixed 
total mass of 2M0. 

where Q{a, A) is the absorption efficiency factor. To calcu- 
late the radial optical depth we use the assumption, that 
the radial extent of the disk is small compared to the dis- 
tance of the inner rim to the star, so that we can consider 
the ray going through the disk material on a constant 
height z. The optical depth on such a ray is given by 



r(z. A) = CdAr = 



3/dArEoQ(a, A) 



4/9dav 27r/io 



exp 



2hl] 



(25) 



In Fig. H we show the optical depth for our standard 
disk as a function of height z/r^ in the disk. Since we keep 
the total disk mass constant, the optical depth is different 
for different grain sizes. What we have plotted is the op- 
tical depth at very short wavelength, using a geometrical 
cross section for the grains. This number is an estimate for 
the optical depth at all wavelengths smaller than 27ra. For 
interstellar grain sizes (0.1 /^m), very high optical depths 
of up to 10^ can be reached, but only at short wavelengths. 
The SED measurement with the longest wavelengths is 
the 850 /im SCUBA measurement, and the 100 /im grains 
show a better estimate for the optical depth at these wave- 
length. Here the standard disk model can reach r = 10'^, 
and with a disk mass of 2Afo, r = lO'' might be feasible. 
As shown in Fig. ^, this optical depth gives a reasonable 
fit to the observed SED even though the flux at 60 pvn. 
and 850 /im is still slightly over predicted. We therefore 
conclude that a massive disk can indeed reach very high 
optical depth. But we also find that the grains should not 
be much bigger than 100 /xm, or the optical depth will 
become too low. 

6.3. Grain settling 

In order to reach the required covering faction, the optical 
depth of the disk must be high not only in the midplane, 
but also at higher z, up to z/r^ « 0.3. Since the disk 
opacity is entirely due to dust grains, but the disk scale 
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Fig. 9. Settling times for dust grains of different size. 




Fig. 10. Ratio of total radiation pressure on a radial col- 
umn of material over the gravitational attraction from the 
binary. Different curves are for different grain sizes with 
a fixed total mass of O.3M0. The dashed line marks the 
blowout limit. 



height is entirely due to gas pressure, gas and dust have 
to be well mixed. From protoplanetary accretion disks it 
is known that dust grains will eventually start to settle 
towards the midplane. The time for a dust grain to settle 
from height z to height zq in a quiescent hydrostatic disk 



is given by Miyake & Nakagawa (1995) 



TT So 1 , z 
2 Pdfl "k 



(26) 



where i^k = \/ GM/r^ is the Kepler frequency in the disk. 

Figure ^ shows the settling times for different grain 
sizes in the standard disk. The timescale varies from 
10^ years for 1 /im grains to 10^ years for 1 cm grains. 
Therefore, at least for grains smaller than about 100 ^va, 
grain settling does not provide a serious constrain to the 
disk model. However, for grains larger than a few cm, the 
lifetime of the disk becomes so short that we would be 
lucky to see HR4049 in this stage at all. 



6.4. Radiation pressure and disk stability 

To estimate at what height the radiation pressure will ac- 
tually be large enough to blow away the entire disk (gas 
and dust), we make the following calculation. We con- 
sider a radial "beam" of material in the disk. The beam 
is all the material outside a surface element AS* in the 
disk wall. The radiation pressure available for that beam 
of material is given by the momentum flux of all the radi- 
ation entering AS* and being absorbed in the beam, thus 
AS'(1 — e"'^) /47rrQC. In order to blow away a part of the 
disk, the radiation pressure must be strong enough to lift 
the entire mass outside of the surface element of the rim. 
This mass is given by p(2:)ArAS'. The ratio of radiative 
and gravitational force is then 



4TTGMcp{z)Ar ' 



(27) 



and (3=1 marks the limit where the radiation force can 
remove the entire beam of material. Inserting the expres- 
sion for the density Eq. (|2^ ) and solving for z under the 
assumption r ^ 1, we find that the critical height Zcrit 
where the disk is blown away by the stellar radiation is 
given by 



-2hn In 



L^, 1 ^/2TTho 



(28) 



Fig. [l^ shows the value of /? for the standard disk and 
3 different grain sizes as a function of height in the disk. 
The positions where the curves cross the line /3 = 1 is the 
location of Zcrit- We can see that for grains larger than 
« 10 fim, the radiation pressure is never large enough to 
remove the disk (gas-|-dust). For low height this is due 
to the large mass in a radial beam of material. For large 
height, the disk becomes optically thin, and an insufficient 
fraction of the radiation is absorbed in the disk. Only if 
the dust in the disk consists mainly of grains smaller than 
10 /im, the parts of the disk with z/vq > 0.8 can be com- 
pletely removed by radiation pressure. 

6.5. Grain drift 

Removing a part of the disk is not the only effect radia- 
tion pressure on dust grains can have. Since the radiation 
pressure only acts on the dust grains, the grains can be 
pushed outwards through the gaseous disk. In parts of the 
disk which are optically thick, all but the innermost grains 
are shielded from direct stellar radiation. The internal ra- 
diation field in the disk is almost perfectly isotropic, and 
the radiation force and corresponding drift motion will 
be small. However, in the optically thin parts of the disk 
and also at the inner rim, grains do receive the full stel- 
lar radiation. In these regions, the velocity at which the 
grains move can be calculated by assuming an equilib- 
rium between the radiation pressure on a dust grain and 
the drag force due to the motion of the grain. In the case 
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Fig. 11. EquiUbrium drift velocity of grains assuming di- 
rect exposure to unattenuated stellar light. Note that the 
drift velocity is given in units tq per year. 

of subsonic drift, the equilibrium velocity is given by (e.g. 



Dominik et al., 1989) 



A^r^c^l^v,^ Eo '''Pbi/^/ ^^^^ 

In Fig. ^ we have plotted the drift velocity as a function 
of height in the disk. In the disk midplane, the drift ve- 
locity is about 10~''ro/yr. Grains can therefore be pushed 
into the disk on a timescale of 10'* years only. This defines 
the tightest constraint on the life time of the disk. Note, 
however, that this limit was calculated under the assump- 
tion of a quiescent disk. Turbulent motions near the in- 
ner disk boundary could be efficient in keeping the grains 
close to the inner boundary. At larger height in the disk, 
the drift becomes very efficient, z/tq — 1/3, grains can 
be pushed through the disk in less than 10^ years. This 
provides an effective truncation mechanism for the disk 
height and may explain the difference between the calcu- 
lated hydrostatic height of the gas disk and the observed 
covering fraction. 

7. Discussion 

The circumstellar environment of HR 4049 is certainly 
unusual in many ways. There are three major constraints 
which are given by the observations, and which must be 
explained by any model: 

1. The UV extinction is highly variable, non-grey and in 
phase with the orbital motion. 

2. The IR excess contains 1/3 of the stellar luminosity. 

3. The spectral shape of the excess is a perfect blackbody 
from 1 /im to 850 /im. 

From (|]) one can immediately conclude that the dis- 
tribution of material is not spherical. Our linc-of-sight 
passes through a region with a strong gradient in extinc- 
tion, probably the upper part of a disk-like structure. The 
grains causing the extinction must be small since large 



grains would imply grey extinction. If a disk-like geometry 
is causing the reprocessing of stellar light, (2) implies that 
the disk must be at least marginally optically thick and 
that it must have a large height. The most difficult part 
turns out to be explaining (3). We have presented two dif- 
ferent scenarios on how a disk may produce a blackbody 
emission: (i) large grains in a marginally optically thick 
disk with very small radial extent and (ii) grains of un- 
known size in a disk with very high optical depth. In this 
discussion we will highlight the virtues and the shortcom- 
ings of both models. 



7.1. Large grains versus small grains 

If the blackbody shape of the IR excess is caused by ther- 
mal emission of large grains (see Sect. these grains can 
only be distributed in a disk. This disk must have a large 
height and it must be confined to at most 10 AU. As dis- 
cussed in Sect. O, it is unlikely that this disk would be a 
gas-poor disk. Moreover, the presence of strong molecular 
bands in the ISO-SWS spectrum of HR 4049 (see Paper 
I) strongly suggests that the circumstellar environment of 
HR 4049 is indeed gas-rich. The gas pressure would then 
support the disk to keep the large scale height. The main 
issue to address is then : where are the small grains ? Some 
small grains are still required for explaining the UV ab- 
sorption as well as the features due to nano-diamonds and 
PAHs in the ISO-SWS spectrum. But there must be an 
efficient mechanism to concentrate the small grains in a 
small region through which our line-of-sight passes, while 
effectively removing them from the regions where most of 
the reprocessing of stellar radiation happens. 

One way to explain this would be to assume that the 
small grains are a completely separate component which 
is not located in the disk, but above the disk or in bipo- 
lar lobes created by a more recent mass-loss event. This 
dichotomy between large and sm all grains is actually con- 
sistent with the conclusions by Joshi et al. ( 1987 ) who 
argued that the polarization they observed points to two 
populations of dust grains, one being very small grains, 
the other being large grains. However, such a geometry 
makes it more difficult to understand the large variations 
of extinction in phase with orbital variations. 

As the gas which is observed through the molecular 
bands is typical for an 0-rich environment (see Paper I), 
it is reasonable to assume that the dust grains which re- 
side in this gaseous reservoir are also 0-rich. If this is true, 
the complete lack of silicate features in the spectrum is an- 
other indication for the absence of small grains in this disk. 
These 0-rich dust grains would then incorporate the re- 
fractory elements lacking in the photosphere of HR 4049. 
If there are mainly large grains in the disk, the dust mass 
as derived from Eq. (|l^) should be a fair estimate for the 
total dust mass. Adopting a standard dust-to-gas ratio 
then yields a total disk mass of a few times 10^^ M0. 
This is somewhat low for the total amount of mass typi- 
cally lost on the AGB. As discussed in this paper, it seems 
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disk erosion due to 
radiation pressure 



Fig. 12. Geometry sketch of the wall model for HR4049. The double star is at the center of the fat disk. The inner 
surface of the disk is heated by the star and produces the blackbody radiation. The material behind the inner surface 
is shielded and cold, including the top, bottom and rear surfaces. The height of the disk is set by strong radiation 
forces which blow away the upper parts of the disk. The direction to the observer is such that we look through the 
edge of the disk. 



unlikely that a large fraction of the dust has escaped the 
system. Moreover, there is no clear-cut evidence (yet) for 
any kind of outflow. The rather low total disk mass would 
then indicate that the binary system prevented HR 4049 
to have a "normal" evolution and leave the AGB prema- 
turely, before the star could experience the most dramatic 
episodes of mass loss. 

Apart from explaining the location of the small grains, 
the main problem with the large-grain model appears to 
be how to support the required disk height for an extended 
period of time. Settling (in a gas-rich disk) and collisions 
(in a gas-poor disk) both limit the lifetime severely. 

The wall model is an appealing alternative to explain 
the blackbody shape of the IR continuum. In this case, 
a standard grain size distribution could be present in 
the circumbinary disk. The blackbody emission would be 
due to the entire surface instead of the individual grains. 
The visible grains near the surface all do have high 
temperatures of about 1200K, and the absence of silicate 
features would be due to the absence of grains in the 
right temperature regime, and to the small temperature 
gradient near the surface. The variations observed in the 
UV absorption - due to small grains - can be tied to 
the vertical density distribution. Indeed, as explained in 



Sect. 6.2, the differential settling times for different grain 
sizes would cause a layered grain size stratification in 
which any large grains would settle towards the midplane, 
while the smallest grains would still be present high 
in the disk, causing the extinction in a beam passing 
through this region. If we assume that the vertical density 
distribution of the small grains is not too affected by 



the settling mechanism, the density distribution of these 
small grains would follow the gas density as given in 
Eq. (p2[). As the small change in viewing angle between 
photometric minimum and maximum causes a large 
change in the amount of absorbed energy, we probably 
view the system with a rather high inclination, skimming 
just past the edge of the disk. 



7.2. The inclination of HR4049 

We can actually get a more quantitative idea about the 
inclination at which we view the system by expressing it 
in terms of the hydrostatic scale height. From Eq. (22) it 
follows that 

dlnp 



din 2: 



''0 



(30) 



To calculate d In z and d In p, we consider the following. At 
photometric minimum, the star is at inferior conjunction 
(i.e. closest to us, poor observers). The line of sight then 
crosses the disk at height i7, thus probing material with a 
density Pmin- At photometric maximum, the star is further 
away from us, and therefore the line of sight crosses the 
disk at + Ai/, this time probing material with a density 
Pmax- From simple geometric considerations it then follows 
that 

(31) 



Rn 



with D the distance between the star at photometric min- 
imum and at photometric maximum and Rq the distance 
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between the center of mass of the binary system and the 
disk. At the same time, the densities are proportional to 
the absorption optical depth so that 



(32) 



and 



d\np dp/p 



Pn 



H 



d\nz dzj Z Pmin 

Finally, Eq. ( ^ ) becomes 



'max ^ \ Ro 




D 
(33) 



(34) 



Using the "grey" optical depths as derived from the UV 



absorption (see Sect. 3.1), i?o=10AU(as in our standard 
mode l in Sect. |) and D --1.18 AU( |yan Winckcl et "al| , 
1995| ) we find that z ~ 2.8hQ. With Hq the pressure scale 



height (Eq. (20)) we then find that the line-of-sight crosses 
the disk at a height H of about 5.3 AU. The inclination i 
is then found by solving 



tan(90° 



H 



Ro - D/2 



(35) 



yielding an inclination of 60°. This shows that the 
variation in the optical depth can be understood in terms 
of the scale height of a hydrostatic disk. It also shows that 
the disk must be truncated to about 2-3 pressure scale 



heights, consistent with what we found in section 6.5 



corresponding to an H/R of 0.58. However, the dynamic 
truncation of the disk means that the density distribution 
in not a perfect Gaussian. The inclination derived above 
is therefore an approximation. 

The Wall model we presented in this paper is there- 
fore not only capable of explaining the blackbody shape 
of the IR continuum, but can also explain in a consistent 
way the variations in the UV absorption tied to the orbital 
variations. In that case, the chemical composition of the 
dust grains and the exact geometric configuration of the 
different grain species observed through either the UV ab- 
sorption or the features in the ISO-SWS spectrum remains 
a matter of debate. The total mass in the circumbinary 
disk which is required to make the optical depths sufh- 



ciently high (see Sect. 6.2) on the other hand might well 
be compatible with the total mass lost during the evolu- 
tion of HR 4049 on the RGB and/or the AGB taking into 
account that no mass has left the system. 
The main difficulty with the Wall model is accurately re- 
producing the exact blackbody shape which we unambigu- 



7.3. The metallicity of HR4049 

We have also derived properties concerning the disk stabil- 
ity. We showed that the central part of a massive disk can 
withstand the radiation pressure acting on the dust grains, 
because the radiation force will only act on the innermost 
grains, while the rest of the disk does not receive direct 
stellar radiation. It is obvious, that a much more thorough 
study of the dynamics of such a disk would be very inter- 
esting. Such a study requires two- fluid hydrodynamics and 
radiative transfer involving large optical depth. We may 
however speculate, that a study of the true dynamics of 
such a disk may provide a natural explanation for the ex- 
tremely low metallicity of IIR4049. Envision the disk filled 
with small grains. The radiation pressure of the star will 
push the grains into the inner rim of the disk, leading to 
very high dust densities there. The dust grains will collide 
with the gas and provide an outward force to the entire 
medium. But since the outer disk does not receive direct 
stellar light, it is not moving radially. Therefore, either a 
large radial pressure gradient or substantial sub-Keplerian 
rotation must balance the radiation force. However, we 



have shown in section 3.5 that the grains will slowly drift 
outwards through the gas, reducing the dust content in the 
innermost layers. Once a layer of gas has become dust free, 
the radiation force acting on this layer will become very 
small. The layer falls inwards, driven either by the gas 
pressure gradient, or by its sub-Keplerian rotation. The 
falling gas will be free of dust and therefore have very low 
metal abundances. If this gas can eventually be reaccreted 
by the primary of HR4049, it will form a metal-poor pho- 
tospheric layer on the star. Since the disk may contain a 
significant fraction of a solar mass, a substantial amount 
of metal free gas is available for this mechanism. 



7.4. Fat circumstellar disks in post- AGB binary stars 

Geometrically thick circumstellar disks are not uncom- 
mon in post- AGB stars. Similar to HR4049 are objects 
with a detected binary with a rather small orbital period 
implying strong binary interaction during evolution. A fa- 
mous example is the Red Rectangle (HD 44179), where 

|1995| ). The Red 



the disk is even resolved (Roddier et al 



Rectangle also shows variability in phase with the or- 
bital motion of the binary, but the variability is entirely 



grey. Waelkens et al. (1996) explain the variability in this 
source with a variable angle under which the stellar light 
is scattered towards the observer. This model also requires 
the presence of a fat disk around the system, similar to 
our proposed wall model. Another example is HD 213985, 
which shows the same photometric behavior as HR 4049 
( Waelkens et al. , 1995| ). The link to the central binary 



characteristics is less clear for objects where t he disc is re- 
ously observe. As was shown in Sect^ O leakage^through ^ol^^^ like the Egg Nebula (|Sahai etg], |l99|). The nature 



the top, bottom and rear surfaces of the disk still produces 
a significant contribution to the flux at long wavelengths, 
even for very high optical depths. 



of the central binary in this object is not yet known. 

However, none of the above sources shows the perfect 
blackbody shape of the excess emission. 



16 



Dominik et al: The dust disk of HR4049 



8. Conclusions 

The Spectral Energy Distribution of HR 4049 shows three 
significant keys which can be used to determine the basic 
properties of its circumstcllar dust: A variable, non-grey 
circumstellar extinction, an infrared excess which amounts 
to 1/3 of the stellar flux and the perfect blackbody shape 
of the excess between 1 /im and 850 /im. We conclude 
that the bulk of the circumstellar matter must be in a 
disk with large scaleheight H/R ~ 1/3. We have con- 
sidered large grains in a marginally optically thick disk 
as a possible explanation and find that such a structure 
is subject to very tight constraints on the radial extent 
and optical depth of the disk, and that the large scale 
height is difficult to support for an extended time. We 
propose, as an alternative solution, the wall model which 
is a fat disk with a radial optical depth in excess of 10^. 
The wall model fulfills many of the requirements, with the 
only significant limitation that the optical depth must be 
very high to suppress excess flux at 850 /im. The wall can 
explain the observed variable extinction in a very consis- 
tent way, while the big-grain model requires a separate 
component for this observation. In the context of binary 
evolution, the formation of a massive circum-binary disk 
can be qualitatively understood as the result of tidal in- 
teraction during the AGB phase of the primary. The wall 
model is therefore an attractive candidate to describe the 
CS environment of HR4049. 

In the near future, the MIDI interferometer at the 
VLT should be able to resolve the disk around HR4049. 
Wavelength dependent interferometric measurements will 
give information about the location of PAH emission and 
continuum, and about the location of the inner disk rim. 
HR4049 has fascinated astronomers now for about 20 
years - there is more to come. 
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